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The Landau— Zener effect is generalized for many-body systems with pairing residual 
interactions. The microscopic equations of motion are obtained and the 14 C decay of 
223 Ra spectroscopic factors are deduced. An asymmetric nuclear shape parametrization 
given by two intersected spheres is used. The single particle level scheme is determined in 
the frame of the superasymmetric two-center shell. The deformation energy is computed 
in the microscopic-macroscopic approximation. The penetrabilities are obtained within 
the WKB approximation. The fine structure of the cluster decay analyzed in the frame 
of this formalism gives a very good agreement with the experimental ratio of partial 
half-lives obtained in special conditions. 

1. Introduction 

Recently, a new method was proposed to approach the fine structure phenomenon 
of cluster- and alpha-decays in odd-nuclei. These decays modes are considered 
as superasymmetric fission processes. Assuming the existence of few collective 
variables, associated to some generalized coordinates describing the nuclear shape 
parametrization, this approach allows to handle approximately the behavior of many 
other intrinsic variables. At any values of the shape generalized coordinates, the 
single particle level scheme and the potential barrier are determined. During the 
whole process and after the disintegration, the single particle occupation proba- 
bilities of the orbitals of interest are computed using the Landau-Zener promotion 
mechanism alone. This method was successfully used to describe the 14 C- and alpha- 
decays 1 ^ 2 ' 3 . Despite the overlook of residual interactions, the agreement obtained 
with data gives a strong experimental support to the formalism. Unfortunately, the 
approach based on the Landau-Zener effect alone is not able to describe the fine 
structure in the case of the disintegration of even-even nuclei. That motivated us 
to develop the model in order to actually obtain the final occupation probabilities 
of the nucleons in different orbitals using the equation of motion derived for pairing 
residual interactions including the Landau-Zener effect. 

The decaying system provides a time dependent single-particle potential in which 
the nucleons move independently. This description is consistent within the essence of 
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the Hartree-Fock approximation where the many body wave function is constrained 
to be Slater determinants and with the Hartree-Fock-Bogoliubov (HFB) problem 
where the function is constrained to be of BCS form all time. In the HFB model, 
the level slippage is treated automatically due to the simultaneous presence of many 
Slater determinants in the wave function. Such calculations are very cumbersome. 
A way to bypass this problem is to use the superfiuid model constructed over a 
single particle potential and introducing the Landau-Zener effect to replace the 
effect given by a part of the residual interactions. 

In the next section, the equations describing the microscopic dynamics of a many 
nucleon system with pairing residual interaction, including the Landau-Zener effect 
applied to the unpaired nucleon, are deduced. This effect allows the transition of 
an unpaired nucleon between levels with the same quantum numbers associated to 
some symmetries of the system. Solving the equations, the response of the system 
is estimated when the nuclear shape parametrization is changed. In the following, 
a nuclear shape parametrization given only by two intersected spheres, the radius 
of the sphere associated to the light fragment kept as constant, will be used. The 
single degree of freedom remains the elongation characterized by R, the distance 
between the centers of the nascent fragments. Sometimes a normalized elongation 
(R — Ri)/(Rf — Ri) will be used for convenience. The parameter R\ being the 
radius of the daughter determined from the condition of volume conservation and 
i?2 = -ilA^ 3 being the radius of the light fragment kept as constant, the next values 
arc determined Ri = R\ — R2 and Rf = Ri + i?2 for the configuration of the initial 
nucleus considered as a sphere and for that of two tangent spheres, respectively. 

2. Formalism 

In order to obtain the equations of motion, we shall start from the variational 
principle taking the Lagrangian as 

C=<4>\H-ih— + H'-\N\<f>> (1) 

and assuming the many-body state formally expanded as a superposition of time 
dependent BCS seniority one diabatic wave functions 

>= e c ^k + n>' w + v ^M a P 1 ° > ( 2 ) 

Each of these seniority one diabatic wave functions being associated to a state 
with an unpaired nucleon in the orbital i, the orbitals j 7^ i being pairwise filled. 
Furthermore 

H(t) = ek(t)(a+a k + a^a- k ) -G^ a+a^a T ai (3) 

fe>0 k,l>0 

is the time dependent many-body Hamiltonian with pairing residual interactions, 

H ' = \ hijitfKuidi - Viaf)(u*af + v*a 3 ) + ( Uj a 3 - Vj af)(u*af + v*aj)} 

(4) 
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is the residual interaction between diabatic levels characterized by the same quan- 
tum numbers associated to some symmetries of the system. This interaction is 
responsible for the existence of avoiding level crossing levels in the case of adiabatic 
states. The presence of this interaction term H' allows the level slippage of an un- 
paired nucleon and allows the possibility that a unpaired nucleon (or hole) located 
on the diabatic level i jumps on the level j due to the interaction hij (t), in the same 
manner as in the classical two level Landau-Zener effect 4 ' 5 . The element hij(t) is 
different from zero only in the regions of the avoided crossings. The term — XN rep- 
resents the well known constraint on the total number of particles. In this picture, 
a~l and a, denote operators for creating and destroying, respectively, a particle in 
the diabatic state i (not, as usual, in the adiabatic state i) and the label % denotes 
the time-reversed orbital conjugate to, and degenerate with, i. Furthermore, e k are 
the single particle energies obtained within the shell model. The blocking correla- 
tions (the variation in A and in the uj 's and vj 's due to the change of the blocking 
level 6 ) are neglected. 

By performing the variation of the Lagrangian in a way similar as in Ref 7 , the 
next system of coupled differential equations are obtained. 

ihpi = ^2p m {KiA* m - K^A m } (5) 

m 

ihkt = 5> m {(2 P; - l)A m - 2 K; (e,(t) - A(t))} (6) 

m 

ifiPm = — y ^ h m j(t)(S m j — Sj m ) (7) 
3 

ihS jm = S jm {-±(\ A m | 2 - | Aj | 2 ) + (e m (i) - €j{t))- 
i[-^+2< + £t-2«;]E I P I A J - 
+ 2Km + % - 2 Kj ] EiPi&t} + 3 E k ?mAhmk(t)S jk - h jk (t)S km ) + \h jm {t)(p 

(8) 

where 

Pm — c rn Crn 
Sjm — Cj Cm 

A m = GJ2k^m K k /n\ 

^m — U l^k^m K k 

K k = U* k V k 
Pk =| V k | 2 

In this description, p. m represents the single particle occupation probability of the 
orbital m, and p k is the pairwise occupation probability of the orbital k. The 
previous system is a generalization of the equations of motions obtained in Rcfs. 
7,s ' 9 and are applied in the case of fine structure of the 14 C emission from 223 Ra. 
This is a rare phenomenon where the theory 10,11 has superseded the experiments 
12 ' 13 . However, the predictions were not able to assess that transitions to the 
first excited state are favored. The system is solved by taking as initial values 
those corresponding to the ground state lin/2 of the parent 223 Ra. The starting 
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elongation is that which gives the minimal value of the microscopic-macroscopic 
deformation energy. 

3. Results Concerning the 14 C Decay of 223 Ra 

The neutron and proton level schemes are obtained with the superasymmetric 
two-center shell model 17 improved in Ref. 1 . The single particle neutron levels are 
displayed in Fig. 1. The 223 Ra has the spin | emerging from lin/ 2 . Adiabatically, 
this unpaired neutron reaches the 2g g / 2 level of the daughter 209 Pb. In the frame 
of this formalism, the fine structure of the 14 C radioactivity can be understood 
by an enhanced transition probability of the unpaired neutron from the adiabatic 
level with f2 = 3/2 emerging from lin/ 2 , to the adiabatic level with the same 
spin projection emerging from 1 ,715/2- The level scheme presented in Fig. 1 shows 
that the lin/2 level reaches adiabatically the 2g g / 2 state, the lj'15/2 reaches the 
1*11/2 state, and the 3^/2 level reaches lj'15/2 state of the daughter 209 Pb. In 
system with cylindrical symmetries, Landau-Zener transitions can be realized only 
between levels with the same spin projection 0. Four avoided level crossings can 




R-RO/(R-Ri 



Fig. 1. Neutron level scheme for 14 C spontaneous emission from 223 Ra with respect the normalized 
elongation. The levels (with $7=3/2) emerging from lin/2, IJ15/2 an d are represented with 

thick lines. 



be observed in Fig. 1 produced at (R - Ri)/{Rf - R t ) ~ 0.75, 1, 1.1 and 1.2. 
Without these interactions, diabatically, the levels lin/2! lj'15/2 an d 3d 5 / 2 attain 1 
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the H11/2, l?'i5/2 an d 2<7g/2 daughter orbitals, respectively. 

In Ref. 8 , a particular form of the equations (5-8) (without Landau-Zener in- 
cluded) were used to study the dumping (or the dissipation) process. The dissipa- 
tion was defined as the flow of energy between collective and intrinsic modes. This 
energy flow must be irreversible. Unfortunately, accordingly with the discussion 
of Ref. 8 , neglecting the Landau-Zener effect, the calculation is microscopically re- 
versible in the sense that if all the coordinates are time-reversed the system retraces 
its path. An unknown fraction of the energy that is identified as dissipated energy 
is in fact collective kinetic energy. Through the presence of the Landau-Zener term 
in the equations of motion the requirement of irreversibility is satisfied. Therefore, 
The Landau-Zener probabilities p m describe a true dissipation process, while the 
variations of pk and values take partially into account the collective kinetic en- 
ergy. In these circumstances, it can be assumed that, after the disintegration, the 
spectroscopic factor can be approximated with only the leading term: 

S m = |< $ | $ m >| 2 OC Pm (10) 

where $ is the BCS wave function given by the equations of motion after the 
scission, $ 0m = a+ Hjj<: m ( u 0j +v 0o ajaf ) | >, with m=2g 9/2 , lin/ 2 , lj'15/2, are 
the orthogonal set of seniority one wave functions constructed on the asymptotic 
overlapped energy levels of the daughter and the emitted nucleus with the ground 
state uoj and vqj values. 

The partial half-life for the state i is obtained with the equation in the WKB 
approximation: 

ft In 2 (2 f Rm I A Y A 2 , \ ft In 2 , 

T "h = [n L V 2 — E * AR ) = SbST"**-' (u) 

where p m is the final occupation probability and has the meaning of a spectroscopic 
factor. E™ is the deformation energy computed in the framework of the Yukawa- 
plus-exponential model 14 extended for binary systems with different charge den- 
sities 15 . This deformation energy is corrected within shell effect computed within 
Strutinsky's method as in Ref. 16 . The deformation energies for the excited states 
are obtained by adding the difference between the single particle energy of the ex- 
cited state and the single particle energy of the nucleon in the ground state. The 
values of K, are 2.10xl0 27 , 1.27xl0 29 and 1.87xl0 31 for transitions to the daugh- 
ter ground state, first excited state and second excited state, respectively. These 
deformation energies are plotted with thick line in Fig. 2. i? g s - is the ground state 
elongation used as starting point of the decay, R l is the channel exit point from the 
barrier and E v is the zero point vibration energy. The effort is mainly focused to 
reproduce the experimental ratio R of the partial half-lives for transitions to the 
first excited state and to ground state, 

R _ T i/2 2 = Pg 9/ 2 cxp(if, 11/2 ) 
7f; 2 /2 Pin/2 exp(X S9/2 ) 

which has approximately the value 0.218. A critical parameter in the model is the 
quantity which characterizes the variations in time of the generalized coordinates. 
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Fig. 2. The potential barrier measured from the ground state of the parent as function of the 
normalized elongation coordinate (R — Ri)/(Rj — Ri). The barriers of the first and second excited 
states obtained by adding the single particle excitation energy are also plotted. 
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In the present work, it is matter of the velocity v of the inter-nuclear distance (or the 
elongation). Some models are able to determine this quantity 18 , but in the present 
work the velocity of the elongation is considered as a fitting parameter and its 
optimal value will be briefly discussed in comparison with previous results. In Fig. 
3, the final single-particle occupation probabilities as function of the internuclear 
distance velocity are plotted. As expected, for large velocities the system becomes 
mainly diabatic and the occupation probabilities of the first daughter excited state 
tends to 1 while that of the ground state reaches 0. For lower velocities, the g.s. 
transitions are favored. In the bottom of Fig. 3, the ratio R is plotted. When 
u« 9 x 10 4 — 3 x 10 5 , R reaches the vicinity of the experimental value. The value 
of v determined in this way is approximately one order of magnitude lower than 
that determined from the tunneling times of Ref. 18 using macroscopic models, 
that means, approximating the effective mass with a value close to the reduced 
mass. However, in the region of avoided crossing levels, the more realistic cranking 
approximation predicts that the inertia increase with orders of magnitude. If the 
effective mass increases, the deformation velocity must decrease in order to conserve 
the available energy. Therefore, the lower values of v are qualitatively consistent 
with the use of the cranking model for the inertia. 

The previous treatment provides a different way to attack the fine structure of 
superasymmetric fission for odd and even nuclei. To our knowledge, it is the unique 
model which can explain the favored transitions to the excited state in the case 
of cluster-decay Also, the inclusion of the Landau-Zencr effect in the equations 
of motion offers a more complete description of dumping phenomena during the 
disintegration. 
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